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ABSTRACT 

We have obtained high-precision interferometric measurements of Vega with the CHARA Array and 
FLUOR beam combiner in the K' band at projected baselines between 103 m and 273 m. The measured 
visibility amplitudes beyond the first lobe are significantly weaker than expected for a slowly rotating 
star characterized by a single effective temperature and surface gravity. Our measurements, when 
compared to synthetic visibilities and synthetic spectrophotometry from a Roche-von Zeipel gravity- 
darkened model atmosphere, provide strong evidence for the model of Vega as a rapidly rotating 
star viewed very nearly pole-on. Our model of Vega's projected surface consists of two-dimensional 
intensity maps constructed from a library of model atmospheres which follow pole-to-equator gradients 
of effective temperature and surface gravity over the rotationally distorted stellar surface. Our best 
fitting model, in good agreement with both our interferometric data and archival spectrophotometric 
data, indicates that Vega is rotating at ~91% of its angular break-up rate with an equatorial velocity 
of 275 km s~^. Together with the measured v sin «, this velocity yields an inclination for the rotation 
axis of 5°. For this model the pole-to-equator effective temperature difference is 2250 K, a value much 
larger than previously derived from spectral line analyses. A polar effective temperature of 10150 
K is derived from a fit to ultraviolet and optical spectrophotometry. The synthetic and observed 
spectral energy distributions are in reasonable agreement longward of 140 nm where they agree to 5% 
or better. Shortward of 140 nm, the model is up to 10 times brighter than observed. The far-UV flux 
discrepancy suggests a breakdown of von Zeipel's Tcgoc g^/^ relation. The derived equatorial Toff of 
7900 K indicates Vega's equatorial atmosphere may be convective and provides a possible explanation 
for the discrepancy. The model has a luminosity of ~37 Lq, a value 35% lower than Vega's apparent 
luminosity based on its bolometric flux and parallax, assuming a slowly rotating star. The model 
luminosity is consistent with the mean absolute magnitude of AOV stars from the VF(H7) — My 
calibration. Our model predicts the spectral energy distribution of Vega as viewed from its equatorial 
plane; a model which may be employed in radiative models for the surrounding debris disk. 

Subject headings: methods: numerical — stars: atmospheres — stars: fundamental parameters (radii, 
temperature) — stars: rotation — stars:individual (Vega) — techniques: interfer- 
ometric 



I. INTRODUCTION 

The bright star Vega (a Lyr, HR 7001, HD 172167, AO V) has been a photometric standard for nearly 150 years. 
iHearnshawl l)1996|) describes Ludwig Seidel's visual comparative photometer measurements, beginning 1857, of 208 
stars reduced to Vega as the primary standard. Tod ay, precise absolu te spectrophotometric observations of Vega 
are available from the far-ultraviolet to the infrared (Bo hlin fc Gi lliland 2004). The first signs that Vega may be 
anomalously luminous appea r in the 1960s after the calibration of the H 7 equivalent width to absolute visual magnitude 
(W^(H7) — Mv) relationship l)Petrielll964D . iMillward fc Walkeil l)1985() confirmed Petrie's findings using better s pectra 
and showed that Vega's My is 0.5 magnitudes brighter than the mean AO V star based on nearby star clusters. IPetrid 
(^964) suggested the anomalous luminosity m ay indicate that Vega is a binary, however the Intensity Interferometer 
measurements bv lHanbu rv Brown et al. (1967) found no evidence for a close, bright co mpanion, a re sult later confirmed 
by speckle observations (McAlister. .198511 . A faint companion cannot be r uled out llAbsil et al J 120 0611. however the 
presence of such a companion would not solve the luminosity discrepancy. iHanburv Brown et alJ 1)196'^ also noted, 
based on their angular diameter measurements, that Vega's radius is 70% larger than that of Sirius. Recent precise 
interferometric measurements show Vega's radius (R = 2.73±0.01 R.^. -Ciardi et al 200L1 to be 60% larger than that 
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of Sirius (R = 1.711±0.013 R©, M =2.12±0.06 M^. lKervella et al.ll2003jl . while the mass-luminosity and mass-radius 
relations for Sirius, Lll.c^ = {MJM^)'^-^^ '^-'^, R/Rq = (M/Mq)0-715±o.035^ yi^y radius for Vega only ~12% larger. 

Since the work of Ivon Zeipell l)1924albD . it has been expected that in order for rapidly rotating stars to achieve 
both hydrostatic and radiative equilibrium, these stars' surfaces will exhibit gravity darkening, a decrease in ef f ective 
temperature from the p ole to t he equator. In the 1960s and 1970s considerable effort (see e.g.. lCollinsl ll963'. '196(tI: 
iHardorp fc Strittmatten.1968 : Maeder fc Pevtremann 197fl; .Collins fc Sonneborn 197 7) was put into the development 
of models for the accurate prediction of colors and spectra from the photospheres of rapidly rotating stars. These 
early models showed that in the special case where one views these stars pole on, they will appear more luminous than 
non-rotating stars, yet have very nearly the same colors and spectrum. The con nection between V ega's anomalous 
properties and the predictions of rapidly rotating model atmospheres was made bv lGravl l|1985lll988j) who noted that 
Vega must be nearly p ole-on and rotating at 90% of its angular breakup rate to account for its excessive apparent 
luminosity. iGravl l|1988f) also noted that Vega's apparent luminosity is inconsistant with its measured Stromgren color 
indices which match that of a dwarf, while the apparent luminosity suggests an evolved subgiant. 

Another anomalous aspect of Vega is the flat-bottom shaped a ppearance of rnany o f its weak metal li nes observed 
at high spectral resolution and very high signal-to-noise (> 2000) f Gulli ver et alJir99lD . The modeling bv lElstd |1992) 
showed that such flat-bottomed or trapazoidal shaped profiles could result from a strong center-to-limb variation in 
the equivalent width of a line coupled with a latitudinal temperature gradient on the surface of the star. Soon after, 
IGulliver et al.l 1)1994(1 modeled these unusual hue profiles together with Vega's spectral energy distribution (SED) and 
found a nearly pole-on («=5.5°), rapidly rotating (T4q — 24 5 km s~^) model to be a good match to these data. 

Since the detection in the infrared of Vega's de bris disk ((Aumann et al.lll984D . much of the attention paid to Vega 
has been focused in this regard (see e.g..l Su et al.12005) . Howe ver, not only has Vega's disk been spatially resolved, so 
too has its photosp here, first by Hanbur y Brown et alJ l|1967|) . though attempts to measure Vega's angular diameter 
go back to Galileo l(Hughesll2001tl . Recent interferometric measurements of Vega show nothin g sig nificantly out of the 
ordinary when compared to a standard models for a slowly-rotating AO V star (Hill et al. 20041 v sin i ~ 21.9±0.2 
km s~^). Specifically, uniform disk fits to data obtained in the first lobe of Vega's visibility curve, from the Mark 
HI int erferometer (M ozurkcwich et al. 2003) at 500 nm and 800 nm and from the Palomar Testbed Interferometer 
(PTI) l|Ciardi et al.ll200ll) in the K-band, show the expected progression due to standard wavelength-dependent limb 
darkening: 3.00±0.05 mas (500 nm), 3.15±0.03 mas (800 nm), 3.24±0.01 (K-band). In addition, the fi rst lobe data 
in th e optical from the Navy Prot otype Optical Interf erometer (NPOI) yield 3.11±0.01 mas (~650nm) l)Ohishi et alJ 
l2004|) . consistent with this picture. iCiardi et al.l l|2001fl note small residuals in their angular diameter fits that may be 
due to Vega's disk. 

Triple amplitude data from NPOI in May 2001 l|Ohishi et al. 2004) sample the second lobe of Vega's visibility curve 
where a gravity-darkening signature should be unambiguous, however these data show the signature of limb darkening 
expected for a non-rotating star, as predicted by ATLAS9 limb-dark ening models (van Hammc 1993|). Most recently, a 
preliminary analysis of second lobe NPOI data from October 2003 ( Peterson et al.ii2004l indicate that Vega is indeed 
strongly gravity darkened, a result inconsistent with Ohishi et al. (2004). Peters on et al. | ( 20061 note that t he NPOI 
Vega data are difficult to analyze due to detector nonlinearities for such a bright star. IPeterson et al.l l(2006}) do see a 
strong interferometric signal for gravity darkening from the rapid rotator Altair with an angular break-up rate 90% of 
critical. Since a sim ilar rotation rate is expected for Vega on the basis of its apparently high luminosity l|Gravlll98§: 
IGulliver et aDIiggl . a strong gravity darkening is expected for Vega as well. 

There is cl early a. nee d for additional high spatial resolution observations of Vega's photosphere to confirm the 
hypothesis of' Gravl l|198 8). confirm the 2 003 NPOI observations, and test the theory of von Zeipel. We have employed 
the long baselines of the CHARA Array l(ten Brummelaar et al.ll2'nn5r on Mount Wilson, together with the capabilties 
of the spatially-filtered Fiber Linked Unit for Optical Recombination (FLUOR,. iCoude du Foresto et al.l l 20 03). as a 
means to probe the second lobe of Vega's visibility c urve at high precision and accuracy in the K-ba nd. O ur Vega 
campaign, part of the commissioning science (McAlist eTet al.i r2005: Mcrand et al. 2005; van Belle et al1 l2006|) for the 
CHARA Array, obtained visibility data on baselines between 103 m and 273 m which clearly show the signature of a 
strongly gravity darkened, pole-on, rapidly rotating star. In this paper we present these data and a detailed modeling 
effort to fit both our inteferometric data and the archival data of Vega's spectral energy distribution. 

We introduce our observations in fj21 Sections SI ^-nd [JHl describe the construction and fitting of one- and two- 
dimensional synthetic brightness distributions to our interferometic data and archival spectrophotometry. A discussion 
of our results follows in We conclude with a summary in S0 

2. OBSERVATIONS 

Our interferometric measurements were obtained using the Center for High Angular Resolution Astronomy (CHARA) 
Array in the infrared K' band (1.94/im to 2.34/im) with FLUOR. Our observations were obtained during 6 nights in the 
late spring of 2005 using four telescope pairs, E2-W2, S1-W2, E2-W1, a nd S1-E2 with maximum basehnes of 156, 211 
251 and 279 m, respectively. The FLUOR Data Reduction Software ()Kervella et al.ll2004t fCoude Du Foresto et alJ 
Il997|) was used to extract the squared modulus of the coherence factor between the two independent telescope apera- 
tures. We obtained 25 calibrated observations of Vega which are summarized in Tabled The (u,w)-plane sampling is 
shown in Figure ^ ^^^^^^^^^ ^_ 

The calibrator stars were chosen from the catalogue of ' Merand et alJ |2005). The CHARA Array's tip-tilt adaptive 
optics system operates at visual wavelengths. Vega is sufficiently bright that it was necessary to reduce the gain on 
the tip-tilt detector system while observing Vega and return the gain to the nominal setting for the fainter calibrator 
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TABLE 1 

CHARA/FLUOR K'-band Vega Measurements 



Projected Position 

Telescope u v Baseline Angle ""Kotal Calibration Star{s) 

No. Julian Date Pair (meters) (meters) (meters) (degrees) X 100 X 100 HD Number 



1 


2453511, 


.261 


E2-W2 


-98 


.941 


23, 


.114 


101, 


.606 


-76, 


.85 


21, 


.1531 


0, 


,8846 


176527 


2 


2453511, 


.313 


E2-W2 


-127, 


.859 


-0, 


.092 


127, 


.859 


89, 


,95 


6, 


,2229 


0, 


,2019 


176527, 


3 


2453511, 


.347 


E2-W2 


-139, 


.876 


-18, 


.250 


141, 


.062 


82, 


,56 


2, 


,6256 


0, 


,0742 


173780 


4 


2453511, 


.374 


E2-W2 


-144, 


.773 


-33, 


.322 


148, 


.558 


77, 


.03 


1, 


.3567 


0, 


,0417 


173780 


5 


2453512, 


.266 


E2-W2 


-103, 


.834 


20, 


.146 


105, 


.770 


-79, 


,02 


18, 


,2301 


0, 


1976 


159501 


6 


2453512, 


.269 


E2-W2 


-106, 


.062 


18, 


.698 


107, 


.698 


-80. 


.00 


16, 


,7627 


0, 


1710 


159501 


7 


2453512, 


.277 


E2-W2 


-110, 


.513 


15, 


.601 


111, 


.609 


-81, 


.96 


14, 


.4223 


0, 


1493 


159501 


8 


2453512, 


.284 


E2-W2 


-114, 


.716 


12, 


.396 


115, 


.384 


-83, 


.83 


12, 


.2229 


0, 


1336 


159501 


9 


2453512, 


.291 


E2-W2 


-118, 


.435 


9, 


.291 


118, 


.799 


-85, 


,51 


10, 


,3873 


0, 


1168 


159501 


10 


2453512, 


.345 


E2-W2 


-140, 


.179 


-18, 


.907 


141, 


,448 


82, 


.31 


2, 


.6399 


0, 


,0741 


173780 


11 


2453512, 


.349 


E2-W2 


-141, 


.068 


-20, 


.951 


142, 


,615 


81, 


.55 


2, 


.3968 


0, 


,0676 


173780 


12 


2453512, 


.356 


E2-W2 


-142, 


.577 


-24, 


.954 


144, 


.744 


80, 


,07 


2, 


,0041 


0, 


,0591 


173780 


13 


2453516, 


.258 


E2-W1 


-141, 


.950 


88, 


.392 


167, 


,221 


-58, 


.08 


0, 


.1040 


0, 


,0059 


159501 


14 


2453516, 


.343 


E2-W1 


-224, 


.986 


25, 


.325 


226, 


.407 


-83, 


,57 


1, 


,2148 


0, 


,0521 


159501, 


15 


2453517, 


.248 


E2-W1 


-132, 


.597 


92, 


.319 


161, 


.569 


-55, 


,15 


0, 


,2426 


0, 


,0194 


159501, 


16 


2453517, 


.288 


E2-W1 


-180, 


.244 


67, 


.502 


192, 


.469 


-69, 


,46 


0, 


,5913 


0, 


,0314 


173780 


17 


2453517, 


.342 


E2-W1 


-225, 


.788 


24, 


.193 


227, 


.080 


-83, 


.88 


1, 


.1066 


0, 


,0670 


173780 


18 


2453519, 


.225 


E2-S1 


-169, 


.006 


-165, 


.745 


236, 


.716 


45, 


.55 


1, 


.1361 


0, 


,0414 


159501 


19 


2453519, 


.252 


E2-S1 


-168, 


.472 


-183, 


.482 


249, 


.095 


42, 


,55 


0, 


,9120 


0, 


,0344 


159501 


20 


2453519, 


.285 


E2-S1 


-161, 


.265 


-205, 


.029 


260, 


.851 


38, 


.18 


0, 


.6047 


0, 


,0259 


159501 


21 


2453519, 


.316 


E2-S1 


-147, 


.913 


-224, 


.292 


268, 


.673 


33, 


.40 


0, 


.5079 


0, 


,0238 


159501 


22 


2453522, 


.270 


E2-S1 


-163, 


.306 


-200, 


.735 


258, 


.773 


39, 


,13 


0, 


,5911 


0, 


,0427 


159501 


23 


2453522, 


.306 


E2-S1 


-148, 


.868 


-223, 


.205 


268, 


.295 


33, 


,70 


0, 


,4518 


0, 


,0241 


159501 


24 


2453522, 


.336 


E2-S1 


-131 


.105 


-239, 


.777 


273, 


.279 


28, 


.67 


0, 


.3788 


0, 


,0199 


159501 


25 


2453538, 


.206 


W2-S1 


56 


.624 


202, 


.948 


210, 


,699 


15, 


,59 


0, 


,9303 


0, 


,0682 


162211 
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Fig. 1. — The sampling of the {u, v)-plane for the CHARA/FLUOR Vega data set. The diamonds represent the monochromatic sampling 
at 2.0 fim. within the K' band. In the K' band, the CHARA baselines E2-W2, E2-W1, E2-S1, and W2-S1 sample the lower first lobe, 
first null, and second lobe of Vega's visibility curve. Two-telescope observations have a 180° ambiguity in the position angle, therefore we 
plot two coordinates, {u,v) and {—u, —v), for each of the 25 data points. These {u,v) points overlay a model for Vega's two-dimensional 
monochromatic Fourier appearance. This squared visibility model is a Fast Fourier Transform (displayed with a logarithmic stretch) of a 
synthetic intensity map of Vega in the plane of the sky (see Figure l3l. 



Aufdenberg et al. 



1.000E~~f 



0.001 - 




100 



150 200 250 

Projected Baseline [meters] 



300 



Fig. 2. — Best fit one-dimensional, symmetric models in comparison with the CHARA/FLUOR data set. The dotted line is a bandwidth- 
smeared uniform disk (xj = 38, Biju = 3.209 ± 0.003 mas) The dashed line is a bandwidth-smeared PHOENIX model atmosphere with 
parameters consistent with a slowly rotating Vega, {T^g = 9550 K, log(g) = 3.95, xi = 20, f^LD = 3.259 ± 0.002 mas), and the solid line 
a bandwidth-smeared analytic limb-darkening model, = (xt = 1.5, (9ld = 3.345 ± 0.006 mas, a = 0.341 ± 0.013). If extended 

emission in the K' band is present at the 1.3% level in the Vega system, these best angular diameters are systematically high by ~3(t (see 
text). 



Stars. Calibrators chosen for this work are all K giants: HD 159501 (Kl III), HD 165683 (KO III), HD 173780 (K3 III), 
HD 176567 (K2 III), and HD 162211 (K2 III). The spectral type difference between the calibrators and Vega does not 
significantly influence the final squared visibility estimate. The interferometric transfer function of CHARA/FLUOR 
was estimated by observing a calibration star before and after each Vega data point. In some cases a different calibrator 
was used on either side of the Vega data point (see Table The inteferometric efficiency of CHARA/FLUOR was 
consistent between all calibrators and stable over each night at ^85%. 

3. ONE-DIMENSIONAL MODEL FITS 

Under the initial assumption that Vega's projected photospheric disk is circularly symmetric in both shape and 
intensity, we have fit three models to the CHARA/FLUOR data set: (1) a uniform disk, where the intensity, assumed 
to be Planckian /(A) = B{Tes = 9550 K, A), is independent of n, the cosine of the angle between the linc-of-sight 
and the surface norma l; (2) an analytic hmb-darkening law, /(/x, A) = B{Tes = 9550 K, A)/x"; and (3) a PHOENIX 
()Hauschildt et al]ll999D mode l radiation field with stell ar parameters {Tcs= 9550 K, log(g) = 3.95) consistent with 
the slowly rotating model thatlBohlin fc Cilli landl l|2004fl show to be a good match to Vega's observed SED. The com- 
putation of the synthetic squared visibilities from these models takes into account the bandwidth smearing introduced 
by the non- monochromatic FLUOR transmission (see § 14.2.11 below) . 

Figure El shows the synthetic squared visibilities from the three models in comparison with the CHA RA/FLUOR 
data. The uniform disk angular diameter we derive is (6'ud = 3.209 ± 0.003 mas) is not consistent with lCiardi et all 
l)2001[l . 9tju = 3.24 ± 0.01 mas. We find this is most likely because we do not assume a flat spectrum across the K' 
band fllter. Regardless, this uniform disk model is poor fit {x1 = 38) because it neglects limb darkening. The limb 
darkening expected for a slowly rotating star should be well predicted by the PHOENIX model, but this model is also a 
poor fit {x1 = 20, 6'ld = 3.259 ± 0.002 mas). The second lobe data indicate Vega is significantly more limb darkened 
compared 1;o this model. The non-physical = ^" model yields a much better fit {xl — 1-5) and a significantly 
larger angular diameter — 3.345 ± 0.006 (a — 0.341 ± 0.013), which suggests the limb-darkening correction in the 
K ' band is ~2.5 tim es larger (4.2% vs. 1.6%) than expected for a slowly-rotating Vega. 

lAbsil et alJ |2006) report that a small fraction, / =1.29±0.16%, of Vega's K' band fiux comes from an extended 
structure, most likely Vega's debris disk. In order to gauge the significance of this extra flux on the photospheric 
parameters derived above, the synthetic squared visibilities are reduced by an amount equal to the square of fraction 



CHARA Array Observations of Vega 



5 



of light coming from the debris disk. At long baselines, the visibility of the debris disk is essentially zero such that: 

n 2 



*^ obs 



(1 — /) ^photosphere + /Vdisk 



(1) 



i0.974K 



photosphere 



The revised fits to ^photosphere are 6Iud 



3.198 ±0.003 [xl = 38) for the uniform disk, 6Ild = 3.247 ±0.002 [xl = 19) 
for the PHOENIX model, and 6Ild = 3.329 ± 0.006 (a = 0.328 ± 0.013, x^=l-4), for the /(/u) = /i" model. The effect of 
removing the extended emission is to reduce the best fit angular diameter for all three models by ~3(t; the correction 
for extended emission is therefore significant. 

4. TWO-DIMENSIONAL MODEL CONSTRUCTION 

In order to physically interpret the strong limb darkening we find for Vega, w e have adapted a computer program 
written by S. Cranmer (private communication, 2002) from lCranmer fc Owockil l)1995j) which computes the effective 
temperature and surface gravity on the surface of a rotationally distorted star, specifically a star with an infinitely 
concentrated central mass under uniform angular rotation, a Roche-von Zeipel model. This azimuthally symmetric 
model is parameterized as a function of the colatitude given the mass, polar radius, luminosity, and fraction of the 
angular break-up rate. 

Each two-dimensional intensity map is characterized by five variables: ^equ, the angular size of the equator, equivalent 
to the angular size as viewed exactly pole-on; uj — il/ricrit, the fraction of the critical angular break-up rate; T^^"^" , 
the effective temperature at the pole; log(.g)poioj the effective surface gravity at the pole; and V'j the position angle of 
the pole on the sky measured east from north. 

Given these input parameters, along with the measured trigonometric parallax TThip = 128.93±0.55 mas 
ijPerrvman et alJll997|) . and the observed projected rotation velocity, v sin i = 21.9±0.2 km s-i l|Hill et all 1200^ . 
the parameterization of the intensity maps begins with 



i?cnu = 107.48- 



(2) 



the equatorial radius in s olar units with both 6'oqu and TThip in milliarcseconds. It follows from a Roche model 
l)Cranmer fc Owockilll995l equation 26) that the corresponding polar radius is 



-^polc — 



3 cos 



TT + cos ^(lj) 

3 



and the stellar mass is 



M 



n 

</polc -Ttpolc 

G 



where G is the universal gravitational constant. The luminosity is then. 



L 



5polc 



(3) 



(4) 



(5) 



where a is the Stefan-Boltzman constant and E is the surface-weighted gravity S ijCranmer fc OwockillT995l equations 
31 and 32), expressed as a power series expansion in w. 



47rGA/ 



1.0 - 0.19696^2 _ 0.094292^4 + 0.33812 



-1.30660W* + 1.8286cjio _ 0.92714^1^ 



(6) 



The ratio of the luminosity to S provides the proportional factor between the effective temperature and gravity for 
von Zeipel's radiative law for all colatitudes -d: 



L 



r^pole 



cff 



5polo 



(7) 



where the gravity darkening parameter, (3, takes the value 0.25 in the purely radiative case (no convection). The 



effective temperature difference between the pole and equator, ATcff, may be expressed in terms of T^^'^'^ and 



AT, 



off 



rynpolc rriCqU 

^eff ~ ^eff 



rxnpolc ( -, 
^Off \ ^ 



27 



7]UJ 



(8) 
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where, 

1] = 3 cos 

The effective gravity as a function of is given by 



TT + cos ^{lj) 



1/2 



—CM 

grW^ijf^ + RWi^sin^f 



(9) 



(10) 



g^{^)^R{i})fl^smi}cosi9 (11) 

where gr and g^) are the radial and colatitudinal components of the gravity field. The colatitudinal dependence of the 
radius is given by 



= cos 

u; sm tf 



n + cos ^ {iv sin ?9) 



(w > 0) 



and angular rotation rate is related to the critical angular rotation rate ^ by 



8 GM 
27 i?^ 



1/2 



At the critical rate (w = 1), i?oqu = l-5i?poio- The inclination follows from 



where the equatorial velocity is 



I = sm 



14. 



V sm z 



-^equ 



(12) 
(13) 

(14) 
(15) 



4.1. Building the Intensity Maps 

For each wavelength, A (185 total wavelength points: 1.9/im to 2.4/im in 0.005/im steps, with additional points 
for H I and He I profiles calculated in non-LTE), an intensity map is computed as follows: Teff{'d) and log{g{'0) are 
evaluated at 90 points from 0° to 90° + i. At each d there are 1024 longitude <p points from 0° to 360° to finely 
sample the perimeter of the nearly pole-on view. For Vega's nearly pole-on orientation, the relatively high resolution 
in ip reduces numerical aliasing when the brightness map is later interpolated onto a uniformly gridded rectangular 
array as described below. 

Each set of spherical coordinates (-R(i9), f) is first transformed to rectangular (x, y, z) coordinates with the 
Interactive Data Language (IDL) routine POLERECSD^. Next, the z-axis of the coordinate system is rotated away from 
the observer by an angle equal to the inclination i (using the IDL routine R0T_3D) and then the x-y plane is rotated 
by an angle equal to V', the position angle (east of north) of the pole on the sky (using the IDL routine ROTATEJCY). 

At each point in the map, the cosine of the angle between the observer's line-of-sight and the local surface normal 

is 

^J■{x,y) = = 

{'^) [sin sin i cos f -I- cos ■& cos i] 



g{'9) 



-gi}{'9) [sin i cos p cos -d — sin -d cos i] 



(16) 



The intensity at each point {x,y) is interpolated from a grid of 170 spherical, hydrostatic PHOENIX (version 13.11.00B) 
stellar atmosphere models (^Hauschildt et al. 1999) spanning 6500 K to 10500 K in Teff and 3.25 to 4.15 in log((7): 

Tj = 6500 -h 250 • j K j = {0, 1, . . . , 16} 
log(5/) = 3.25 + 0.1 / = {0,1,...,9}. 

Four radiation fields, /(A, [i) evaluated at 64 angles by PHOENIX, are selected from the model grid to bracket the local 
effective temperature and gravity values on the grid square, 

T, <T,B{'d)<Tj+i 
gi < gi'd) < gi+i. 



There is a typographical error in equation (5) of lColhn^ 119631) which is not in the paper's erratum lCoUinslll964ll : 



GM 
He 



should 



be = ^^r, where luc the critical angular rate, and Re is the equatorial radius at the critical rate. 

^ The coordinate transformation routines used here are from the JHU/APL/SIR IDL library of the Space Oceanography Group of the 
Applied Physics Laboratory of The Johns Hopkins University. 
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The intensity vectors I\{fJ,) are linearly interpolated (in the log) at fj,{x,y) around the grid square, 

lT^hiT„gi,fi{x,y)) 

Ix' ^ h{Tj, gi+i, ii{x,y)). 

Next, the intensity is bilinearly interpolated at the local Toff and \og{g) for each (a;, y) position in the map: 

Ix{x,y)=Ix[Tcs{x,y), g{x,y), fi{x,y)] 

= (l-a)(l-fe)/r + a(l-fe)/r 

+a6/ii + (1 -a)6/f (17) 

where 



a=(T^ff{x,y) -Tj)/{Tj 



+1 



T, 



b^{g{x,y) -gi){gi+i ~ gi) 

Finally, a Delaunay triangulation is computed (using the IDL routine TRIGRID) to regrid the intensity map I\{x, y), 
originally gridded in 'd and ip^ onto a regular 512x512 grid of points in x and y. The coordinates x and y have the 
units of milliarcseconds and correspond to offsets in right ascension and declination on the sky (Aa,A5) relative to 
the origin, the subsolar point. 

4.2. Synthetic Squared Visibility Computation 

Due to the lack of symmetry in the synthetic intensity maps, we evaluate a set of discrete 2-D Fourier transforms 
in order to generate a set of synthetic squared visibilities comparable to the CHARA/FLUOR observations. The first 
step is to compute the discrete Fourier transform for each wavelength at each of the spatial frequency coordinates (u, v) 
corresponding to the projected baseline and orientation of each data point (see Table^. The mean (u,w) coordinates 
for each data point, in units of meters, are converted to the corresponding spatial frequency coordinates (uk,Vk) in 
units of cycles per arcsecond for each wavelength A^. The Fourier transform. 



V^{u,v) 



5A/A(x,2/)e^2^("^+"^) dxdy 



— oc ^ — oc 



(18) 



is approximated by the integration rule of Gaussian quadrature fe.g.. IStroud fc SecrestillQB^ iPress et alJll99^ 



V^{uk,Vk)' 



' N N 

^ ^ AjSkhixi, yj) cos{2n{ukXi + Vkyj)) 
.i=i j=i 

' N N 

Ai ^ AjSkhixt, yj) siii{2Tr{ukXi + Vky^)) 



(19) 



where 5*^ is the wavelength discretized value of the instrument sensitivity curve 5*^, and Ai, Aj and Xi, yj are the 
weights and nodes of the quadrature, respectively. For our square grid, the x— and y— coordinate nodes and weights 
are indentical. The 2-D Gaussian quadrature is performed with a version of the IDL routine INT_2D modified to use an 
arbitrarily high number of nodes. The intensity at wavelength fc, Ik{x,y), is interpolated at {xi,yj) from the regular 
512 X 512 spacing to the quadrature node points using the IDL routine INTERPOLATE which uses a cubic convolution 
interpolation method employing 16 neighboring points. The synthetic squared visibility is normalized to unity at zero 
spatial frequency by: 

N N 



^fc(0,0) 



'^A,J2^jSkIk{xi,yj) 

,i=i j=i 



(20) 



We find N = 512 provides the degree of numerical accuracy sufficient in the case of a 2-D uniform disk (right circular 
cylinder) to yield values in agreement with the analytic result. 



Vk{uk,Vk) 



2M7r9^ul + vl)/{ne^ul + vl) 



(21) 



(where Ji is the first order Bessel function of the first kind, 9 is the angular diameter of the uniform disk and B is 
the projected baseline), to better than 1% for > 10^^. We use the IDL function BESELJ for our Ji computations. 
For V'^ < 10~^, near the monochromatic first and second zeros, the numerical accuracy of the quadrature deteriorates 
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Fig. 3. — A synthetic brightness map (hnear stretch) of Vega for our best fitting parameters: lo = 0.91, 9cqu = 3.329 mas, T^'^^ = 10150 
K, log((;)poie = 4.10. For this model, Vega's pole is inclined 5° toward a position angle of 40° and 'x' marks the subsolar point. The labeled 
intensity contours are relative to the maximum intensity in the map. 



to 10% or worse. The bandwidth-smeared V'^ minimum is at ~ 10~'^, so we find this quadrature method yields 
squared visibihties which are s ufficiently accurate for our task, however observations with an even larger dynamic 
ranee ( Perrin & Ridgwavl l2005|) will require more accurate methods. 

To test the 2-D Gaussian quadrature method in the case where no analytic solution is available, we computed the 
2-D Fast Fourier Transform (IDL routine FFT) of a brightness map (see Figure I^J. First, we compared the results of 
the 2-D FFT to the analytic uniform disk, equation 1)21(1 . To reduce aliasing we find it necessary to "zero pad" the 
brightness map. With 12-to-l zero padding (the 512 x 512 brightness map placed at the center of a larger 6144 x 6144 
array of zeros) we find the 2-D FFT has very similar accuracy to the 512-point Gaussian quadrature: better than 1% 
down to > 10^^ inside the second null. For the brightness map shown in Figure O the 2-D FFT and Gaussian 
quadrature methods agree to better than 0.5% down to V'^ > 10^^, inside the second null. We find the computational 
time required to evaluate equation (|19|l at 25 (wfc, Vk) points for 185 wavelengths is 6 times faster than the evaluation 
of the 185 zero-padded 2-D FFTs. 

4.2.1. Bandwidth Smearing 

Once we have computed V^{uk,Vk) for the 185 wavelength points, we proceed to compute the bandwidth-smeared 
average squared visibility V{B,Xq)^, 

This integral is performed by the IDL routine INT_TABULATED, a 5-point Newton-Cotes formula. The term is 
included so that the integral is equivalent to an integral over wavenumber (frequency) where 

-1 _ /r 

is the mean wavenumber. This simulates the data collection and fringe processing algorithm used by FLUOR. In 
the discretized integrand, V{B,Xk)^ is equivalent to V^{uk,Vk) where B — 206264. SA^^u^ + v^, with A^ in units of 
meters and Uk and Vk in units of cycles per arcsecond. 

4.3. Synthetic Spectral Energy Distribution Construction 

To construct synthetic SEDs for Vega from the Roche- von Zeipel model, 170 radiation fields were computed from 
the same model grid used to construct the K' band intensity maps. The wavelength resolution is 0.05 nm from 100 nm 
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to 400 nm and 0.2 nm from 400 nm to 3 /zm and 10 nm from 3 fiia to 50 /im. The higher resolution in the ultraviolet is 
needed to sample the strong line blanketing in this spectral region. From the resulting grid of radiation fields, intensity 
maps are computed fsee 14. 1|) and the flux is computed from 



This 2-D integral is performed with the IDL routine INT_TABULATED_2D (version 1.6) which first constructs a Delaunay 
triangulation of points in the ■dip-plane. For each triangle in the convex hull (defined as the smallest convex polygon 
completely enclosing the points), the volume of the triangular cylinder formed by six points (the triangle in the plane 
and three points above with heights equal to the integrand) is computed and summed. For computing the flux from 
the intensity maps, a coarser sampling in and ip (20x40), relative to that needed for the visibility computations, 
is sufficient for better than 1% flux accuracy. The numerical accuracy was checked by computing the SED of a 
non-rotating star (w = 0) and comparing this to a single effective temperature SED from a 1-D atmosphere. The 
interpolation and integration errors result in a flux deficit of less than 0.7% at all wavelengths relative to the 1-D 
model atmosphere. 



The computation of each intensity map, the Fourier transforms, and the bandwidth-smearing for each set of input 
parameters {9equ, w, T^g^'^, log((7)poic, V') is too computationally expensive to compute synthetic squared visibilities 
many hundreds of times as part of a gradient-search method over the vertices of a 5-dimensional hypercube. Therefore, 
we must proceed with targeted trial searches to establish the sensitivity of xt to each parameter after first establishing 
a reasonable range of values for each parameter. 

The parameter 0cqu is a physical angular diameter related to a uniform disk fit by a scale factor depending on the 
degree of gravity and limb darkening, which in turn depends on the parameters w, log((7)poio, and T^^^'^, in order of 
importance. As shown above, a limb-darkening correction of 4% is significantly larger than the ^1.5% value expected 
for a normal AO V star at 2/im (Davis ct al. 2000). The analytic limb-darkening model fit is sufficiently good that we 
take 9cqu — 3.36 mas as a starting value. This corresponds to ii!equ — 2.77 R© from equation ||2J). 

Our starting value for lo is based on the assumption that Vega's true luminosity should be similar to that slowly 
rotating AO V stars. Vega has an apparent luminosity, assuming a single effective temperature from all viewing angles, 
of 57 Lq based on its bolometric flux and the parallax. In the pole-on rapidly rotating case, we would see Vega in 
its brightest projection. According to Millward & Walker (1985) the mean absolute visual magnitude. My, is 1.0 for 
spectral type AO V. With a bolometric correction of —0.3, this translates to i = 37.7 L0. From equations Q and © 
we expect uj > 0.8 in order to account for the luminosity discrepancy assuming a minim um p olar effective temperature 
of 9550 K, based on the non- rotating model fits to Vega's SED (Bohhn & GiUilanJ 1200^ . The range of effective 
temperatures and surface gravities for the model atmosphere grid described in ii4.1l sets our upper rotation limit at 
ui < 0.96. For lu > 0.8, ATcs > 1300 K (see equation |51l , thus T^^^'^ must be greater than 9550 K to compensate for 
the pole-to-equator temperature gradient and to reproduce the observed SED. So, given a mean apparent T^s of 9550 
K, a rough estimate of T^^^^ is 9550 K -I- iATcs = 10200 K. We therefore limit the polar effective temperature to the 

range 10050 K < T^g" < 10350 K. 

The relationship between oj and the true luminosity, through equations (0,(111), and is independent of the polar 
surface gravity, yet we can constrain log((7)poic by assuming Vega follows the mass-luminosity relation we derive for 
the slowly rotating Sirius, L/Lq = (M/M©)^ '^^''-^. Here we assume Vega's rapid rotation has no significant effect on 
its interior in relation to the luminosity from nuclear reactions in its core. Assuming L ^ 37.7 Lq from above, the 
mass-luminosity relation yields M=2.3±0.1 Mq. As lo increases, i?poie decreases relative to i?cqu, therefore choosing 
M = 2.2 Mq and uj — 0.8 provides a lower limit of log(g)poie = 4.0. For lower polar gravities, Vega's mass will be 
significantly lower than expected based on its luminosity, nevertheless we choose a range log((7)poic values from 3.6 and 
4.3 in order to check the effect of the gravity on our synthetic visibilities and SEDs. 

Lastly, the position angle of Vega's pole, ■0, should be important if Vega's inclination is sufficiently high and its 
rotation sufficiently rapid to produce an elliptical projection of the rotationally distorted photosphere on plane of the 
sky. Previous measurements l)Ohishi et alJl200'^ J(^jjj(^^i_j1^jjj|2p( 11[l fin d no evidence for ellipticity. Preliminary results 
from the NPOI three-telescope observations of iPeterson et al.l l|2004ll suggest an asymmetric brightness distribution 
with ^ = 281°. 



For the grid search we compute the reduced chi-square xt for a set of models defined by 6'equ, w, T^g log((7)poio, 
and Ip, adjusting O^qu slightly (<0.3%) to minimize xl for each model (see below). Figure 0] shows a xl niap in the 
Lo — plane for a range of 6loqu values with T^"^" = 10250 K. log(g)poie = 4.1. We find a best fit of xl = 1-31 for 
parameters uj — 0.91, ^cqu = 3.329 mas, and -0= 40°. A direct comparison of this model with the squared visibility 
data is shown in Figure 

The F test provides a la lower limit on uj at ~0.89. For uj < 0.89, the synthetic V'^ values are generally too high 
across the second lobe because the model is not sufficiently darkened towards the limb. Correspondingly, the upper 




(24) 



5. TWO-DIMENSIONAL MODEL FITTING 



5.1. Initial Parameter Constraints 



5.2. CHARA/FLUOR Data: Parameter Grid Search 



10 



Aufdenberg et al. 



150 



100 




tn 

CD 
CD 
l_ 

CP 

CD 

o 

CL 

M— 

O 

Q> 
CD 
C 
< 

c 

'in 
o 
Q. 



0.86 0.88 0.90 0.92 0.94 
Fraction of Critical Angular Rotation Rate, cj 

Fig. 4. — A contour plot of Xv i'l oj — 1/1 plane for r^'° = 10250 K and log(g)poi(, = 4.1. The labeled contours denote the lower and 
upper IfT range, and a 2a contour, from the F test. The 'x' marks the best fit, x^=1.31, while the brightest region has a x2 = 3.25 (see 
Figure ISl. 
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Fig. 5. — a) The CHARA/FLUOR data (error bars) plotted as a function of projected baseline (for a range of azimuths, see Table [Tl 
together with the best fitting Roche-von Zeipel synthetic squared visibilities. Model parameters: <.j=0.91, 9cqu = 3.329 mas, Tj^'^ = 10250 
K, log(g)poie = 4.10. The best fit xi = 1-31. b) Deviations of the best-fit model from observed squared visibilities. The dotted and dashed 
lines indicated the Icr and 2a deviations. 



limits on lu are constrained because the synthetic values are generally too low across the second lobe, due to very 
strong darkening toward the limb for uj > 0.93. In addition, the upper limit on w is a function of ip because the 
projected stellar disk appears sufRciently more elliptical, even at low inclinations j ~ 5°, as the model star rotates 
faster. The data from the nearly orthogonal E2-W1 and E2-S1 baselines constrain models with lu > 0.92 to limited 
range of position angles, but these data provide no constraint on -0 at lower u values where the star is less distorted. 



CHARA Array Observations of Vega 



11 



TSS = 10250 K 

iog(g)poie = 4.1 



_L 




_L 



3.40 



3.38 



3.36 



0) 

■+-' 

V 

E 
g 

Q 



0.84 0.86 0.88 0.90 0.92 0.94 0.96 
Fraction of Critical Angular Rotation Rate, w 



3.34 



3.32 



3.30 



TS5 = 10250 

iog(g)poie = 4 



b) 



_L 



_L 



0.84 0.86 0.88 0.90 0.92 0.94 0.96 
Fraction of Critical Angular Rotation Rate, 



Fig. 6. — Constraints on model parameters from the CHARA/FLUOR data, a) The reduced chi-square values xf, from the Roche-von 
Zeipel model fit to the squared visibility data as a function the fraction of the critical angular break-up rate, lu = 0/Ocrit> for fixed values 
of the polar effective temperature T^g^" and polar surface gravity log(g)poie. The dashed line denotes the la confidence region for uj from 
the F test for 24 degrees of freedom relative to the best fit at xi = 1.31. For each ui, xi values are plotted for 18 position angles ip (0° 
to 170° in 10° steps, see Figure b) The relationship between the best fit equatorial angular diameter dcqu at each u! for the range of 
position angles. The dashed lines provide an estimate for the la range in u! and the corresponding range in the equatorial angular diameter. 



-^Gqu/-^pole ^ 1.24. 

As CO increases so does the darkening of the Unib due to the increasing larger pole-to-equator effective temperature 
difference. As a result, the best fit 0equ value increases with ui because the effective "limb-darkening" correction 
increases. The best fit values for 6'oqu and w are therefore correlated. To establish this correlation, we estimated the 
best fitting 0equ value for a given ui without recomputing the brightness map and Fourier components. While each 
intensity map is constructed for a fixed 0cqu value, we can approximate the squared visibilities for models with slightly 
(< 0.5%) larger or smaller ^equ values as follows. A small adjustment to due to a small adjustment in 6'oqu, assuming 
the physical model for the star is not significantly changed and the model changes relatively slowly with wavelength, 
is equivalent to computing at a larger (smaller) wavelength for a larger (smaller) value of ^cqu- So, for a given 
projected baseline, we linearly interpolate (in the log) V^{u, v) at A = ^k{Sfit/9equ), a wavelength shift of 10 nm or less. 
Near the bandpass edges, the instrument transmission drops to zero so there is no concern about interpolating outside 
of the wavelength grid with this scheme. The normalization, equation (|20|l . must be scaled by the (0fit/6'cqu)^ to 
compensate for the revised surface area of the star. After one iteration, setting ^cqu — ^fitj recomputing the Fourier 
map and refitting the data, the best fit ^oqu value is within 0.25% of that found with the estimated model values. 

Figure EJl shows the xi values from Figure 0] projected on the w axis, with a spread of values for the 18 position 
angles at each uj value. This shows again that for the range 0.89 < uj < 0.92 there is no constraint on the position 
angle of the pole. The corresponding best fit 0oqu values are shown in Figure ^jp. The equatorial angular diameter 
is constrained to the range 3.32 mas < 6'oqu < 3.34 mas. The best fit to the CHARA/FLUOR data is insensitive 
to T'g^g^''. This is because ATeff, which determines the overall darkening, is quite sensitive to cu, but not T^^^° (see 
equation |(SJ)). Thus, we cannot usefully constrain T^g^'^ or ip from the CHARA/FLUOR data. As for the surface 
gravity, varying log(f;)poic over what we consider the most probable range, 4.1±0.1, does not significantly effect the xt 
minimum. Models with log((7)poic values from 3.9 to 4.3 all fall within Icr of the best fit. The best fit O^qu values are 
essentially independent of T^g^'^ between 9800 K and 10450 K and weakly dependent on log((7)poio between 3.8 and 
4.3; all best fit ^equ values fall well within the la range established in Figure El 
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Fig. 7.— A contour plot of xl for the SED fits in the oj - Tfg''' plane. The uj range is limited to the la range from CHARA/FLUOR 
fits (see Figure|SJ. The polar surface gravity is fixed at log(g)pojc = 4.1. The labeled contours denote the Icr and 2a regions from the F 
test. The 'x' marks the location of the best fit model, x^=8.7. 



Here we compare our synthetic SEDs to the absolute spectroph otometry of Vega. Specifically, we compare our models 
to the data-model composite SED^" of l)Bohlin fc G illiland''2004^ which includes International Ultraviolet Explorer data 
from 125.5 to 167.5 nm, HST Space Telescope Imaging Spectrograph data from 167.5 to 420 nm, and a specifically 
constructed Kurucz model shortward of lUE and longward 420 nm to match and replace data corrupted by CCD 
fringing in this wavelength region. To facilitate this comparison, first the synthetic spectra were convolved to the 
spectral resolution of the observations (A/ A A ~ 500) and then both the data and convolved synthetic spectra were 
binned: 2.0 nm wide bins in the UV (127.5 nm to 327.5 nm, 101 bins) and 2.0 nm bins in the optical and near-IR 
(330.0 nm to 10080 nm, 340 bins) for a total of 441 spectral bins. 

Figure 13 shows the xt in the uj — T^^^^ plane. These two parameters, apart from the angular diameter, most 
sensitively affect the fit to the observed SED. There is a clear positive correlation between lu and T^g^°. This makes 
sense if one considers that a more rapidly rotating star will be more gravity darkened and require a hotter pole to 
compensate for a cooler equator in order to match the same SED. Following this correlation, it is expected that a 
continuum of models from {oj — 0.89, T^f^''^ = 10150 K) to {oj = 0, T^g^° = 9550 K) will provide a reasonable fit to the 
SED since the non- rotating ATLAS 12 model of Kurucz fits the observed SED quite well (Bohlin & Gilhland 2004). 
However, wc did not consider models with u; < 0.88 in the SED analysis because such models are a poor match to the 
CHARA/FLUOR squared visibility data set as shown above. In other words, although the ATLAS 12 model provides 
a good fit to the observed SED, it fails to predict the correct center-to-limb darkening for Vega. 

The best fit synthetic spectrum is shown in Figure |H1 Considering the complexity of this synthetic SED relative to 
a single T^s model, there is generally good agreement (±5%) between our best fit model and the data longward of 300 
nm, apart from larger mismatches at the Paschcn and Balmer edges and in the Balmcr lines. Longward of 140 nm, the 
model agrees with the observations to within ±10%. At wavelengths below 140 nm, as measured by the lUE, the data 
are up to a factor of 2 lower than predicted. Our best fit yields x^=8.7. The overprediction below 140 nm has only a 
small effect on the synthetic integrated flux between 127.5 nm and 10080 nm, 2.79xl0~^ erg cm~^ s~^, which is within 
1.2ct of the value derived from an integration of the observed SED, (2.748±0.036) x 10~^ erg cm"^ s~^. The equatorial 
angular diameter derived from this SED fit, ^oqu = 3.407 mas, differs from the best fit to the CHARA/FLUOR data, 
^oqu = 3.329 mas, by 2.4%, a value within the uncertainty of the absolute flux calibration. 

6. DISCUSSION 

The best fit stellar parameters, based on the model fits to the CHARA/FLUOR data and archival spectrophotometry 
in !J3 are summarized in Table [21 As discussed in Sj^l the effect of extended K' band emission in the Vega system, if 

^'^ ftp://ftp.stsci.edu/cdbs/cdbs2/calspec/alpha_lyrjtis_002.fits 
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a) A comparison between the SED of IBohlin fc Gillilandl 120041) and our best fitting (xt = ^■'^) rapidly rotating SED model 
= 0.91, T^g^'^ = 10150 K and log((ji)poic = 4.10. The differences between this model and the data in the b) region at shorter 
wavelengths observed by the lUE and the c) region observed by the HST Space Telescope Imaging Spectrograph at longer wavelengths. 
At w avenumbers 1/A < 2.38 ^m~^ the "observed" SED is represented by a closely fitting Kurucz model spectrum f see IBohlin fc GillilanJ 
200l. 



TABLE 2 

Fundamental Stellar Parameters for Vega 



Parameter 


Symbol 


Value 


Fraction of the angular break-up rate 




0.91 ±0.03 


Equatorial angular diameter (mas) 


^equ 


3.33 ±0.01 


Parallax (mas) 


""hip 


128.93 ±0.55 


Equatorial radius (R0) 


^cqu 


2.78 ±0.02 


Polar radius (R0) 


^polc 


2.26 ±0.07 


Pole-to-equator T^g difference (K) 




2250+"*°° 


Polar effective temperature (K) 


-'cff 


10150 ± 100 


Luminosity (Lq) 


L 


37 ± 3 


Mass (Mq) 


M 


2.3 ±0.2 


Polar surface gravity (cm s~^) 


log(9)pole 


4.1 ±0.1 


Equatorial rotation velocity (km s"-^) 


Vcqu 


270 ± 15 


Projected rotation velocity (km s"^) 


V sin i 


21.9 ±0.2 


Inclination of rotation axis (degrees) 


i 


4.7 ±0.3 



Reference 



CHARA/FLUOR fit 
CHARA/FLUOR fit 
Perrvman et^^ ^^2£) 
Equation 121 
Equation |3J 
Equation jSJ 

Fit to spectrophotometry ^ohlii^^^illiland^224) 
Equation j^J 

(L/Lq) = (Af/M0)4-27±O.2O girius) 
Equation jl} 
Equations 1131 and 1151 
Hill et a l. (2004) 
Equation 1141 



unaccounted for, is to increase the apparent angular diameter of Vega slightly, by ~0.3%. Correcting for this effect 
via equation the best fit equatorial diameter is shifted systematically lower by 0.3% (0.01 mas) to the range 
3.31 mas < 6'cqu < 3.33 mas. We find all other parameters in Table [51 are uneffected by the extended emission within 
the error bars given. The best fit range for the fraction of the angular break-up rate, 0.89 < oj < 0.92, sensitive to the 
amplitude of the second lobe, is unaffected by the extended emission because the V"^ correction is quite small there, 
AV"^ < 0.0003, relative to the first lobe where the correction is up to 20 times larger. 

One parameter which stands out is our large pole-to-equator effective tem perature difference, ATeg = 2250^30 ° K, 
relative to previous spectroscopic and spectrophotometric studies of Vega l)Gulliver et alJlT99^ iHill et al.ll2004tl for 
which AToff falls into the range 300 to 400 K. Our larger AToff yields a much cooler e quatorial effective temperature, 
T^^"" = 7900t^!^EJ K, than most recently reported for Vega, 9330 K fH ill et alJl20()l . The amplitude of the second 
lobe visibility measurements as observed by CHARA/FLUOR is well fit by strong darkening toward the limb. In the 
context of the Roche- von Zeipel model, such darkening requires a large pole-to-equator T^s gradient. Consequently, 
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Wavenumber [/^m"''] 



Fig. 9. — Left: A comparison between the SED from lBohlin & : Gilliland I2004J) {lUE and i75T observations supplemented by a slowly 
rotating model spectrum both below 127.5 nm and longward of 420 nm) and two rapidly rotating models for Vega's SED, one viewed from 
an inclina tion of 5° (nearly pole on) and one viewed from an inclination of 90° (equator on), from an integration of two intensity maps via 
equation 1241 for these inclinations. Right: A comparison of the best-fit brightness distributions for Vega with inclinations of 5° (top) and 
90° (bottom). For the equator-on view, the poles appear 10% fainter than the pole-on view due to limb darkening. 



TABLE 3 
A Model Equatorial Photospheric 
Spectral Energy Distribution for Vega 
FROM 1020.5 A to 40AtM (R=500). 



Wavelength 




Flux (Fx) ^ 


(A) 




(erg cm~^ s"-' A"'-) 


1.020500000000000EH 


-03 


1.68027E-(-03 


1.021000000000000EH 


-03 


1.65680E-I-03 


1.021500000000000EH 


-03 


1.62296E-I-03 


1.022000000000000EH 


-03 


1.57629E-I-03 


1.022500000000000EH 


-03 


1.51370E-I-03 



''The complete version of this table is in the elec- 
tronic edition of the Journal. The printed edition 
contains only a sample. 

''The flux at a distance d from Vega in the equa- 
torial plane is the flux from column (2) multiplied 
by the ratio (Rcqn/d)-^ , the ratio squared of Vega's 
equatorial radius to the distance, or (2.78/d)^ 
when d has the units of solar radii. 



we predict that Vega's equator-on SED (that is viewed as if i = 90° and integrated over the visible stellar disk, see 
equation H24I) ') has a significantly lower color temperature and overall lower flux, particularly in the mid- ultraviolet 
where the flux is lower by a factor of 5, as shown in Figure El A debris disk, aligned with Vega's equatorial plane 
as suggested by our nearly pole-on model for the star and the recent observations of a circular disk in the mid-IR 
ijSu et alJl2005() . should see a significantly less luminous, cooler SED than we see from the Earth. In the literature to 
date, modeling of the heating, scattering, and emission of Vega's dusty debris disk has assumed an irradiating SED 
equal to the pole-on view of Vega (see e.g.,|Absil et al. 2006; Su et al. 2005). Our synthetic photospheric equatorial 
spectrum for Vega is tabulated in Table|3| It should be interesting to investigate how our predicted equatorial spectrum 
used in such modeling will affect conclusions regarding the amount of dust and the grain-size distribution in the debris 
disk. 

Several of Vega's fundamental stellar parameters (ATcfr, Vcqu, i) we derive differ significantly from those derived by 
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IGuUiver at al.1 l)1994|) and iHill at al.1 l|2004|) from high-dispersion spectroscopy. Regarding ATcff, both spectroscopic 

studies find w ~ 0.5, while we find w = 0.91±0.03. These two w values along with the corresponding T^g^° values, 9680 
K and 10150 K, in equation Q yield ATcff values of 350 K and 2250 K. The reason the uj values differ is at least partly 
linked to inconsistent parameters used in the spectroscopic studies. As noted in HilLet al. (2004), the Gulliver et a ^ 
l) 1994ft study finds a low value for the polar gravity, log((ji)poio — 3.75, which yields a mass for Vega of only 1.3 4 Mq 
and an inclination inconsistent with the expected equatorial velocity. The equatorial velocity of IHill et al.l l)2004l) . Voqu 
= 160 km s~^, is not consistent with their other parameters {uj = 0.47, log(g)poic = 4.0, i?cqu = 2.73 R0, 7.9°) which 
should yield instead Vcqu — 113 km s~^ and j=ll.l°. Values Vcqu — 160 km s^^ and i — 7.9° are recovered if a; = 
0.65, which corresponds to X^qu/K^rit = 0.47. It is possible to confuse uj with Kqu/Krit- The two are not equivalent: 

TT + cos^"^(ti;)l 



For u! — 0.65, the corresponding ATcff — 757 K, not 350 K. Therefore, there appears to be a mismatch between the 
Vcqu and ATcff values used in the most recent spectral analyses and this suggests the spectral data must be reanalyzed 
with a consistent model. A. Gulliver (private communication, 2006) confirms that Hill et al. (2004) did confuse uj with 
Kiqu/V^rit and this group is now reanalyzing Vega's high dispersion spectrum. Our best fit value for uj, derived from the 
interferometric data is appealing because, together with our derived polar effective temperature, it yields a luminosity 
consistent with that of slowly rotating AO V stars. A more slowly rotating model for Vega will have a warmer equator 
and an overall higher true luminosity too large for its mass. Therefore, it seems that less rapidly rotating models for 
Vega do not offer an explanation for the apparent over-luminosity with respect to its spectral type. 

Our best fit model, while it provides self-consistent parameters within the Roche-von Zeipel context, has several 
discrepancies, most notably producing too much flux below 140 nm relative to the observed SE P. The limitation s of the 
LTE metal line blanketing for modeling Vega in the ultraviolet have recently been explored bv lGarcia-Gil et all ([2005). 
They find that in the UV the line opacity is generally systematically too large in LTE because the over-ionization 
in non-LTE is neglected. Our best model flux below 140 nm is already too large, so a fully non-LTE treatment is 
not expected to improve this discrepancy. The Wien tail of Vega's SED will be the most sensitive to the warmest 
colatitudes near the pole. In our strictly radiative von Zeipel model, SEDs with T^"^" < 10050 K produce too much 
flux in the optical and near-IR, so simply lowering T^^f^'*' will not solve the problem, the temperature gradient must 
differ from the Teff oc g^^^ relation. The equatorial effective temperature we derive, 7900 K, may indicate that Vega's 
equatorial region is convective. If so, von Zeipel's purely radiative gravity darkening exponent, /3 = 0.25, will not be 
valid near the equator. A more complex model, where the gravity darkening transitions from purely radiative near 
the pole to partially convective near the equator, may be the next approach to take. Such a temperature profile may 

allow for a cooler T^g^'^, reducing the fiux discrepancy below 140 nm, while still matching the observed optical and 
near-IR fluxes. Such a gradient must also improve the match to the Balmer and Paschen edges and the Balmer lines. 

7. SUMMARY 

We have demonstrated that a Roche-von Zeipel model atmosphere rotating at 91±3% of the angular break-up rate 
provides a very good match to K' band long-baseline interferometric observations of Vega. These observations sample 
the second lobe of Vega's visibility curve and indicate a limb-darkening correction 2.5 times larger than expected for 
a slowly rotating AO V star. In the context of the purely radiative von Zeipel gravity darkening model, the second 
lobe visibility measurements imply a ^^22% reduction in the effective temperature from pole to equator. The model 
predicts an equatorial velocity of 270±15 km s~^, which toge ther with th e measured usin i yields an inclination of 
i ~ 5°, conflrming the pole-on model for Vega suggested bv iGravl ()1988D to explain Vega's anomalous luminosity. 
Our model predicts a true lum inosity for Vega of 3 7zb3 Lq, consistent with the mean luminosity of AO V stars from 
W(Rj) - Mv calibration (Millw ard fc Walkerl[T985j) . We predict that Vega's spectral energy distribution viewed from 
its equatorial plane is significantly cooler than viewed from its pole. This equatorial spectrum may significantly impact 
conclusions derived from models for Vega's debris disk which have employed Vega's observed polar- view spectral energy 
distribution, rather than the equatorial one, which seems more appropriate given our observations. 
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